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Abstract 

The Large Hadron Collider is shown to have great scope for a light charged Higgs discovery, in 
the context of the CP-violating type-I two Higgs doublet model. This scenario with similar masses 
of if ± and W was suggested by the puzzling departure from charged current lepton universality 
found in the LEP data. With the lightest neutral Higgs mass set to 125 GeV, the charged-neutral 
Higgs associated production mechanism can cause a significant excess in the rubb events over a 
vast range of tan/3 as long as the Higgs mixing pattern avoids a few limiting cases. Thanks to 
the low mass, the charged Higgs loop can play a striking role in neutral Higgs decays into 77, 
thereby compensating for a suppressed gluon-gluon fusion rate. The effect of scalar-pseudo-scalar 
mixing on loop-induced Higgs signals is also discussed. 
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I. INTRODUCTION 



There are many motivations to introduce an extended Higgs sector. Examples are, includ- 
ing obsolete ones, fine-tuning in the weak-scale Higgs mass [l[ suggesting supersymmetry 
[2|, CP- violation jl, |iL_grand unification (HI, extended gauge symmetry jgj], strong CP prob- 
lem Q, string theory vacuum stability (jf, light fermion mass/weak-scale hierarchy (To| . 
fermion flavour structure 11], neutrino mass 12], dark matter (see e.g. fl3l . l3|). cosmic rays 
lH , baryogenesis [l6[ , inflatio n [l7Tj , novel collider phenomenolog y [l8| , as well as anomalies 
seen in precision (see e.g. fl9l . |20| |) and accelerator (see e.g. (2lL |22j) experiments. In this 
article, we entertain the possibility that an intriguing outcome from the LEP experiment is 
in fact pointing to an extended Higgs sector, which may be directly inspected at the Large 
Hadron Collider (LHC). 

Specifically, we take as our framework the two Higgs doublet model (2HDM) of type-I 
with Mh± in the vicinity of the H^-boson mass. The two main points of this paper are the 
following. (1) The charged Higgs produced in association with the lightest neutral Higgs 
can be discovered as an excess in the pp — > rvbb process at the LHC. Moreover, this could 
happen even at the early stage with y/s = 7-8 TeV in a substantial portion of the parameter 
space. (2) With the lightest neutral Higgs mass assumed to be ~ 125 GeV, the charged 
Higgs one-loop effect can enhance B(Hi — > 77) enough to produce the excess observed at 
the recent LHC analyses j23-26[. This enhancement is particularly welcome since the Hi 
production at hadron colliders is generically suppressed due to Higgs mixing. 

This light charged Higgs scenario was previously advocated in j2l[, whose contents we 
recap here. The LEP Elecroweak Working Group has performed an analysis of W decay 
branching ratios. In their report, B(W —> rv) appears to be higher than B(W — > ev) 
and B(W —> fiu), to the level of 2.8 standard deviations j27|. Should this difference be 
real, it would violate the charged current lepton universality, that is predicted by the gauge 
invariance of the Standard Model (SM) and has been confirmed in many indirect tests. 

As a resolution of this puzzle, one of the authors pointed out that a significant portion of 
the apparent lepton non-universality can be attributed toproduction of light charged Higgs 
pairs which predominantly decay to the rv final states [21] . One of the most natural ways 
to realise a H ± that has mass ~ mjy and is thus light enough to be on-shell produced at 
LEP, is to employ the type-I two Higgs doublet model (2HDMI). In this model, all the 
if ± -fermion couplings are suppressed like 1/ tan/3, the ratio of the two vacuum expectation 
values (VEVs). With the choice of a high enough tan/3, this enables H ± to escape from 
the severe flavour changing neutral current (FCNC) constraints as well as any process that 
probes the if ^fermion interaction such as Z — > bb, t — > bH ± , and B tv. 

Presuming the above LEP anomaly to be real, the following points make this scenario 
more appealing. First, invoking a light charge Higgs (discussed also in j28[) is the only 
proposed solution that does not spoil the other precision lepton universality tests using 
lepton or meson decays j29|. Second, the requirement of a light charged Higgs plus b — > 57 
singles out type-I out of the four types of 2HDMs 3(| in which tree-level FCNC is forbidden 
by a Z 2 symmetry Note that these assumptions already fix many aspects of the model 
that otherwise offer more than one option, i.e. the Yukawa coupling structure, the charged 
Higgs mass, and the viable range of tan/3. These combine to predict a distinct set of collider 
physics signatures. 

In this paper, we shall consider CP-violation in the Higgs sector, which can in general 
occur in a class of 2HDMs. In this way, we can explore Higgs mixing patterns which received 
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relatively less attention. A practical advantage is that one can analyse both the scalar and 
the pseudo-scalar Higgses in a unified manner: switching between the two sectors reduces 
to tuning the Higgs mixing angles, which determine the composition of each neutral mass 
eigenstate. 

Within the context of the Minimal Supersymmetric Standard Model (MSSM), there are 
earlier studies of CP- violating Higgs phenomenology. In particular, the general CP- violating 



potential and mass spectrum can be found in the paper of A. Pilaftsis and C. Wagner [31 



where the importance of the scalar-pseudo-scalar transitions in Higgs phenomenology was 
emphasised. In the same article, the couplings of the charged Higgs boson to CP-violating 
Higgs bosons were also discussed. 

This article is organised as follows. We begin by giving an overview of the model in 
the next section. Section II I II is a summary of our search strategy, which leads to the main 
phenomenological results presented in section (TV] Section [V] is devoted to the diphoton 
signal of the lightest neutral Higgs. We deliver our conclusion in section IVIl In addition, we 
present: the ghost Lagrangian for our setup of 2HDM in appendix [AJ the relation between 
our parametrisation and perturbativity in appendix |Bl and differential cross-sections of our 
main process in appendix O 



II. SETUP 



Model 



Let us give a brief description of the type-I 2HDM that we have chosen. The type of a 
2HDM refers to one of the multiple ways to organise the Yukawa couplings. There are four 
types in which quarks or leptons of each charge couple to only one Higg s doublet 3pj . so 
that neutral Higgses are naturally prevented from mediating FCNC |32j. To this end, we 
impose the softly broken Z 2 symmetry, under which the two Higgs doublets $i )2 transform 
as 

($ 1 ,$ 2 )^(-$ 1 ,+$ 2 ), (1) 

whereas the fermion fields remain unchanged. This allows only the following Yukawa cou- 
plings, 

-C Y = L L <f> 2 Y e e R + Q L $ 2 Y d d R + Q L e $* 2 Y u u R + h.c, (2) 



between $ 2 and fermions, so all the fermions acquire mass from the VEV of $2 [ID 
The scalar potential including the soft Z 2 -breaking terms reads 



Ai 



A, 



:$t$,) 2 + A 3 ($ t 1 $i)($^2; 



+ A 4 ($I$ 2 )($2$ 



■2^2 
1 

2 L 



A 5 ($I$ 2 ) 2 + h.c. 



(3) 



+ m2 2 



(*5*2)}- 



- - {tt4 [m^l^Hh.c. 

In unitary gauge, the Higgs doublet fields are expanded around the minimum like 

i sinBH + 



73 ^ 



—i cos/3 H + 
-^{v 2 + r] 2 + i cos/3 T] 3 ) 



(4) 
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One can take a basis of $12 such that v\ and v 2 are both real. Their quadrature sum is 
subject to the constraint, vf + v 2 = v 2 , where v = 2~ 1 / 4 G^ 2 is the SM Higgs VEV. Their 
ratio is parametrised by 

tan/3 = v 2 /vi, (5) 

in terms of the pseudo-scalar and charged Higgs mixing angle /3. 

We introduce CP- violation in the Higgs sector by assuming ImA 5 7^ and lmm 2 2 ^ 0, 
which enables the neutral scalar and pseudo-scalar components to mix together. The real 
symmetric mass-squared matrix Ai 2 with respect to the basis (771,772,773), is diagonalised by 
the mass eigenstates 

H i = R ij rj j , i,j = 1,2,3, (6) 

like 

RM 2 R r = dmg(M 2 1 ,Ml,M 2 ), (7) 

using a 3 x 3 orthogonal matrix R. It is parametrised in terms of the mixing angles 0:1,2,3 i n 
the form, 

Ci C 2 Si c 2 s 2 

-(ciS 2 S 3 + S 1 C 3 ) CiC 3 -SiS 2 S 3 C 2 S 3 
— C\ S 2 C3 + S1 S 3 — (Ci S 3 + Si s 2 c 3 ) c 2 c 3 

where s, = sinaj and q = cosaj. The neutral Higgs mass eigenvalues are assumed to be in 
the order, 

M 2 < M 2 < Ml (9) 



R 



(8) 



This determines the following physical domain of the mixing angles [33 



- tt/2 < a 1)2)3 < tt/2. (10) 

The two CP- violating mass matrix elements are related by [34| 

(M 2 ) 13 = tan/3(7W 2 ) 23 , (11) 

which translates into the constraint, 

M 2 = M 2 R 13 (R 12 tan/3 - R n ) + M 2 R 23 (R 22 tan/3 - R 21 ) 
3 R 33 (R 31 -R 32 tan/3) ' 1 j 

We shall remain in the parameter volume that is compatible with the mass ordering ([9j) in 
conjunction with (112j) . in order to avoid double counting of physically identical parameter 
sets. 

Two more dimensionful quantities can be derived from the scalar potential: the charged 
Higgs mass, 

M| ± =/i 2 -^ 2 (A 4 + ReA 5 ), (13) 

and the auxiliary parameter, 

(j, = Rem 2 2 /sin2/3, (14) 

that sets the mass scale of charged/neutral Higgs particles except Hi. 

We have the eight physical quantities determined from the Higgs potential, 

Mi,M 2 ,M H ±,ai,a 2 ,a 3 ,n,ta.n/3. (15) 
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In this work, we opt to invert this dependency and express the parameters appearing in ([3]) 
in terms of those in (fl5j) 34J. The mass parameters are then given by 

m u = Aii>i + (A 3 + A 4 + ReA 5 — 2u) v 2) (16a) 
m 2 22 = \ 2V l + (A 3 + A 4 + ReA 5 - 2v) v 2 } (16b) 
Imm] 2 = lm\ 5 viv 2 , (16c) 

with v = fi 2 /v 2 and RemJ 2 fixed by (1141) . The quartic couplings are 

A i = ^t c i c 2^i + (c lS2 s 3 + s lC3 ) 2 M 2 2 

c j v (17a) 

+ (cis 2 c 3 - S!S 3 ) 2 M 3 2 - sfjfj 2 }, 



1 



A 2 = ^[s^M 2 + (cic 3 - sis 2 s 3 ) 2 M 2 



V (17b) 
+ (cis 3 + sis 2 c 3 ) 2 M 3 2 - c 2 ^ 2 ], 

A 3 = -^^{c lSl [c 2 M 2 + (s 2 s 2 - c 2 )M 2 

+ (s 2 c 2 - sl)Ml] + S2 c 3 s 3 (c 2 - s\){Mt - M 2 2 )} (17c) 

+ 1(2M 2 ± -//), 
v 

A 4 = ^[^M 2 + c\s\Ml + c 2 c 2 M 3 2 + /i 2 - 2M 2 ± ], (17d) 
ReA 5 = ^{-s\Ml - c\s\Ml - c\c\Ml + /x 2 ), (17e) 



ImA 5 = ^{cp[c 1 c 2 s 2 M l - c 2 s 3 (cis 2 s 3 + s 1 c 3 )M 2 

+ c 2 c 3 (sis 3 - cis 2 c 3 )M|] + s /3 [sic 2 s 2 M 1 2 

+ C 2 S 3 (CiC3-SiS 2 S 3 )M 2 -C 2 C 3 (CiS 3 + SiS 2 C 3 )M 3 2 ]}, 



(17f) 



where s@ — sin/3, eg = cos/3. 

In this approach, the solutions for Ai,...^ may turn out to be non-perturbatively large, 
unless one takes care to set the input parameters in ffT5|) to sensible values. To avoid entering 
a nonsensical regime, we shall check that the sizes of the Higgs quartic couplings meet the 
conditions, 

|Ai|<47r, z = l,..., 5. (18) 

Another theoretical constraint on the quartic couplings that is closely related to pertur- 
bativity arises from the tree-level unitarity in Higgs-Higgs scattering. We also take this into 
account using the formulation in 35] . As we find in the numerical analysis, this perturbative 
unitarity is always satisfied as long as we require ffTB"]) . within the parameter space explored 
in this work. 

The interactions of Higgs bosons with other particles are described by the Lagrangian, 

2 

AC = ^(D^tfiD^i) + £ Y , (19) 
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HiUj Uj 


-i(m Uj /v)(R i2 / sin/3 - iR i3 cot/37 5 ) 


Hidjdj 


-i(m dj /v)(R i2 / sin/3 + iR i3 cot/37 5 ) 




V2(cot(3/v)V jk (m Uj P L - m dk P R ) 




-i(mi. /v) (R i2 1 sin/3 + iR i3 cot/37 5 ) 




-V2 (cot/3 /v)m h P R 


H t W+W~ 


i(2rny V /v)(Rn cos/3 + R i2 sin/3)^ 



HiZ^Zy i(2m 2 z /v)(Ra cosfi + R i2 sin/3)^ 

HjH+W- i(m w /v)[(Rii sin/3 - R i2 cos/3)i - R i3 ]{p + p% 

TABLE I. Feynman rules. The momentum directions are indicated in figure [21 The chirality 
projection matrices are defined by Pl/r = (1 T 7 5 )/2- 

with £y from (j2J). One expands the covariant derivatives, replaces $12 by (j4]), and diago- 
nalises the fermion and Higgs mass terms, to arrive at the Feynman rules in table HI They 
will be used to calculate and interpret the results that we present in the following sections. 

A remark is in order regarding the parametrisation in the CP-conserving limit. One can 
enforce CP- conservation simply by setting a 2 = a 3 = 0. This eliminates mixing between the 
CP-odd and the CP-even degrees-of-freedom, as can be seen in (JSJ). The mixing between 771 
and r] 2 with the angle a = ct\ — tt/2, results in the lighter and the heavier CP-even neutral 
eigenstates, h and H. The CP-odd neutral Higgs A is identified with 773. In this case, the 
relation (llip becomes trivial and therefore it no longer constrains Mf to be determined as a 
function of the other parameters. As a result, we are left with two less mixing angles and one 
more free mass parameter in the CP-conserving Higgs sector than in the CP- violating case. 
Conventionally, these seven free parameters are chosen to be: M h , M H , Ma, M H ±,a, (3, \i 
0. 



B. Parameter space 

Having established the framework, one should fix the range of parameters to study. For 
future references, we first collect the values which were chosen for the reasons explained 
below: 

M H ± = 86 GeV, M x = 125 GeV, M 2 = 200 GeV, \i = 100 GeV, tan/3 = 5. (20) 

This set shall be used by default throughout the analysis unless specified otherwise, in 
combination with one of the two Higgs mixing configurations, 

. . . . ((—0.6,0.1,0.5), benchmark point Pi, . . 

(sincci, sma 2 , sma 3 ) = < (21) 

l( 0.0,0.8,0.5), benchmark point P 2 . 

The charged Higgs mass is more or less predetermined by the LEP data. The lower 
bound is provided by the direct search, and the upper limit comes from the requirement 
that the H ± pair production should restore lepton universality in leptonic W decays. We 
take M H ± which reduces the apparent discrepancy between B(W — > rv) and the average of 
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B{W — > eu) and B(W — > pv) by 1 a [2l|. The latest update from the OPAL collaboration 
reports the direct search limit, M#± > 82 GeV, for Ma> Mh± at 95% confidence level (CL) 



(37J. In view of the recent evidences from ATLAS [25] and CMS |26|, we set the lightest 
neutral Higgs mass to 125 GeV. The heavier Higgs mass parameter is taken to be fairly 
higher than M\. Non-degenerate mass spectrum is needed for CP- violating Higgs mixing. 
The last dimensionful parameter p is set to a representative value of the weak scale order. 

The chosen M 2 is a mass at which the SM Higgs production is strongly constrained at the 
LHC Not far from there lies also M3 < 300 GeV. In this range, the rates of the WW and 
ZZ channels are limited to be roughly below half their SM predictions (38[. In our model, 
the upper bounds must be reinterpreted as arising from processes mediated by H2 and H% 
whose interactions are different from those of the SM Higgs. For this, we require that 

a(gg -»■ H 2 , 3 WW/ZZ) < 0.5 a{gg -»■ h SM -»■ WW/ZZ). (22) 



A potentially important constraint is provided by the p parameter [39(. With the Higgs 
mass spectrum assumed in this work, one can check that the one-loop correction from the 
Higgs bosons to p is small enough to pass the electroweak precision tests [2l|, 

S3- 



For safety of a light charged Higgs, it is crucial to choose a hig h enough tan/3. The 
strongest bound is given by b — > s~f, which requires tan/3 > 4 [2l|. This range of tan/3 
also suppresses Bit — >■ bH + ) below 4%. This maximal top-quark branching fraction coin- 
cides roughly with the numerical values of the 95% CL limits from the light charged Higgs 
searches at the LHC j4l[. However, these LHC constraints on t -> bH + are based on 
the assumption that B(H ± — > tv) = 1, and therefore become weaker in our setup where 
B(H ± tv) = 0.7. In principle, one could use any tan/3 larger than the FCNC bound, as 
long as it lets the charged Higgs decay inside the detector. In practice, it is bounded from 
above by perturbativity through f JT7|) . Obviously, this upper limit is a function of the other 
parameters. In appendix (Bj it is shown that p can always be adjusted so that perturbativity 
holds. One should not regard this as a fine-tuning but rather an artifact of the employed 
parametrisation. If one started by setting the parameters appearing in and then derived 
the quantities in (I15p . then one could trivially satisfy perturbativity by keeping A^...^ inside 
the range (11 8p . Even in this reverse approach, one could naturally obtain a high tan/3, as 
one can see from the example of inert doublet model where an exact Z2 symmetry leads to 
infinite tan/3 (l3l |. 

Among the three remaining Higgs mixing angles, we scan ot\ and a.2 which determine the 
composition of the lightest neutral Higgs. In (JS}, one can notice that sina^ is the fraction 
of the CP-odd component in H\. We fix sina 3 to 0.5 to have an intermediate level of 
CP- violation. For illustrative purposes, we select two benchmark points, labelled Pi and P2 
in (12 ip . The fraction of A in Hi is chosen to be small at Pi and large at P 2 . 



III. LIGHT H± SEARCH STRATEGY 

Given the model and the input parameters, we consider which charged Higgs production 
mechanism is most useful. We first calculate the cross-sections of the following standard 
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channels: 



qq -»■ W ->• #1^, (23a) 

gg j j + j j + WH*, (23b) 

^-►W + iZi-^W + Wtf*, (23c) 

qq^Z + Hi^Z + WH±, (23d) 

^ iZ< Wtf*, (23e) 

gb-*tH~, (23f) 

gg,qq^tb + H + . (23g) 

The last two entries are meant to imply their conjugate processes as well. The final state 



in Q23g[ ) can be made through either top pair production or gluon-gluon fusion [42 



In figures [H, we plot the cross-sections against tan/3 for the two benchmark points. One 
observes the general tendency that each production rate is depleted for high tan/3. As one 
can expect from table HJ processes (123fjl and Q23gQ die out as tan/3 grows, while the others 
approach their individual asymptotic cross-sections even though many of them are small. 

An outstanding exception is the production associated with a neutral Higgs (123ajl . which 
is not suppressed even for high tan/3. The diagrams for this process are shown in figure [2l 
which scale according to the H i -H ± -W vertices in table [B In the Hi-H ± -W interaction, 
the pseudo-scalar coupling is proportional to sina 2 and independent of /3, and the scalar 
coupling is proportional to cos «i cos 0:2 i n the limit of tan/3 — > 00. Combining these two 
behaviours, one can expect that the HiH ± production becomes efficient unless |sinai| is 
large and | sin a 2 1 is small. 

The results plotted in figures [1] do not depend on \i. Therefore, they remain valid even if 
one adjusts \i for the perturbativity of A^...^ in the high tan/3 regime. This means that the 
production mechanism fl23al) is of the most general interest among the displayed channels, 
and we shall focus on it in what follows. Note that the same type of diagram was previously 
considered in the contexts of the MSSM j43| and a fermiophobic Higgs scenario jiij . 

The next step should be to select the decay products of H ± and H±. The two dominant 
branching ratios of H + are 

B(H + -> t + v t ) = 0.71, B(H + -> cs) = 0.27, (24) 

as long as M H ± < 135 GeV, beyond which the charged Higgs decays mediated by a virtual 



top becomes non-negligible [45j. These branching fractions are essentially independent of 
any of the input parameters. This is a notable feature of the 2HDMI, stemming from the 
universal scaling behaviour of the if ± -fermion couplings: they are all proportional to cot/3, 
as shown in table [B As we are mainly interested in a charged Higgs that is light enough 
to have been produced at LEP energies, we can regard (I24I) as good approximations and 
simply choose the decay product with the highest rate, rv. 

As for Hi, we plot its branching fractions in the vicinity of M\ = 125 GeV, in figures [31 
One finds that the bb mode is dominant at both benchmark points. This tends to be the case 
unless both sinai and sin 02 are vanishingly small. At the point where sinai = sin«2 = 0, 
the tree-level ifi-fermion coupling vanishes, causing its fermiophobia. Considering these 
dominant decay modes of H ± and Hi, we are going to concentrate on the final state, rvbb 
for numerical simulation. 

One can also notice that the branching fractions of the 77 and gg modes can differ 
significantly from their individual SM values. Albeit not directly related to the charged 
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FIG. 1. Cross-sections of the single charged Higgs production mechanisms plotted against tan/3 for 
the selected benchmark points and centre-of-mass energies. The curves are independent of \x. 



Higgs search, these changes greatly affect the phenomenology of H\. We shall come back to 
this issue in section [V] 
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FIG. 2. Feynman graph for the charged Higgs production in association with a neutral Higgs. 

IV. LHC PROSPECTS OF if± SEARCH 

The main process of our concern is pp — > W — > H ± Hi, followed by the decays, H — > tv 
and Hi — > bb at the LHC. Among these three subprocesses, the charged Higgs production 
and decay revealed their outstanding stability under variation of tan/3 over a wide range, in 
the preceding section. 

We shall envisage two experimental conditions at the LHC, which we call the "early" 
and the "full" searches. The "early search" denotes the centre-of-mass energy of the current 
run, y/s = 7 TeV, and the integrated luminosity of L = 10 fb -1 . We assume that the "full 
search" will deliver L = 100 fb" 1 at y/s = 14 TeV. 

The kinematics of the rvbb state can be fully reconstructed with high efficiency as it 
has only one neutrino. We assume a r-reconstruction efficiency ~ 100% and a 6-tagging 
efficiency ~ 70%, plus the following set of selection cuts: 



which we choose considering the operation of ATLAS. Assuming these cuts, we calculated 
the SM background^] The cross-section is a(pp — > t ± z/66) B g — 55 fb at yfs = 7 TeV and 
~ 110 fb at v^= 14 TeV. 

Along with the associated production of H Hi, we also study a similar process with H ± 
replaced by W, i.e. pp — > W — >■ WH 1 . As we shall see, the parameter space covered by this 
Higgs-strahlung is complementary to the region accessible to the H ± Hi state. We select 
the final state, l^vbb with I = e, /i, and assume a jy-reconstruction efficiency ~ 100%. The 
cuts in (1251) are used again, and the same cut on r is applied on the transverse momenta of 
the light leptons, Pri^) > 25 GeV. Using these parameters, we obtain the Higgs- less SM 
background cross-section to be a(pp — > / ± ^66)bg — 70 fb at yfs = 7 TeV and ~ 140 fb at 



p T (b),p T (T ± )>25 GeV, 
\rj(b)\, 1^)1 < 2.5, 



pr{v) > 20 GeV, 

\V(V)\ <4-5, 
M(b,b) = Mi ± 10 GeV, 



(25) 



v/i = 14 TeV. 



1 Including contributions from a SM Higgs with mass of 125 GeV. 
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FIG. 3. Branching ratios of the (lightest) neutral Higgs plotted against its mass. The SM case is 
shown in panel [3aJ Panels I3bl and [3c] are respectively for the benchmark points Pi and P2 in the 
2HDMI. The parameters are set as in (|20p unless stated otherwise. 



We used LanHEP 46 



parton distribution functions [48 



to implement the model, and CalcHEP [47[ plus the CTEQ6M 
to perform the phenomenological analysis. The Passarino- 



Veltman functions appearing in the effective vertices of Hi-g-g and ifj-7-7 are calculated 
by calling the LoopTools package (49| . 

We first check the rates of our main processes with sin 02 = sin 03 = 0. In this CP- 
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FIG. 4. Production cross-section of each neutral Higgs in association with a charged Higgs via 
a virtual VF-boson in the CP-conserving case. The masses are Mh = Ma = 125 GeV, Mh = 
200 GeV. For each channel, the lighter and the darker curves represent ^fs = 7 TeV, 14 TeV, 
respectively. In panel Hbl the hatched region is excluded by LEP at 95% CL [37| . and each curve 
is for «i — f3 that maximises the cross-section. 



conserving case, one has the lighter and the heavier CP-even Higgses, h = Hi and H = H 2 , 
and the CP-odd Higgs, A = H 3 . The production cross-sections of h, H, and A, each 
in association with a charged Higgs, are shown in figure |4a] as functions of a\ — fi. For 
each type of neutral Higgs, there are three curves that respectively represent H~, H + , 
and H ± productions in the order of increasing thickness. The higher H + rate relative to 
H~ reflects the fact that the LHC is a pp collider, not pp. The shape of each curve is 
essentially determined by the dependence of the corresponding W-H ± -H i coupling on a\ 
and j3, presented in table [D These patterns will be useful for understanding the simulation 
results to follow. From the plot, one can expect O(10 4 ) events that occur through a charged- 
neutral Higgs associated production with the "full luminosity" at 14 TeV, unless the mixing 
angle is such that the given channel is highly suppressed. 

Our analysis is of relevance to a light charged Higgs, even if one is not interested in a 
scenario where charged Higgs bosons have been on-shell produced at LEP energies. In this 
regard, it is legitimate to consider changes that would be caused by varying M H ±. Obviously, 
the above production rates will be reduced as M H ± grows, as is shown in figure [4bJ We 
take the range of M H ± up to ~ 135 GeV, beyond which the changed Higgs decay pattern 



undergoes a qualitative change (45J, thereby invalidating our assumptions given in (1241) . For 
M H ± ~ 135 GeV, each cross-section decreases roughly down to half the value in panel HU 
These curves illustrate the phase-space suppression that takes place for higher M H ±. We 
emphasise that this reduction of statistics is the only essential information that one needs to 
understand how the significances would change in the following quantitative analysis with 
Mh± = 86 GeV, if one took a different Mn±. 
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(a) V« = 7 TeV, L = 1Q fb _1 (b) -Js = 14 TeV, L = 100 ffcT 1 

FIG. 5. Significance contours for the processes, pp — > H Hi — > r^vbb (solid lines) and 
pp — > W Hi — > l^ubb (dashed lines), on the (sin ai, sin a 2 ) plane in the CP-violating case. The 
parameters are set as in (|20p unless stated otherwise. Light grey (green online), dark grey (red 
online), and black correspond to the significances £ = 2,3,5, respectively. The medium-shadowed 
regions (red online) are excluded by theoretical constraints. The dark-shadowed region (blue on- 
line) is ruled out by current limits on -#2,3- The light-shadowed region (green online) allows for 
discovery of both processes. 



For the simulation of charged Higgs search, we turn on CP-violation allowing for mixing 
among all the three neutral Higgses, h, H, and A. In figures |5l the contours of the signif- 
icance S are shown for each process, which have been evaluated by means of the method 



described in [50|. There are two complications in comparison to the CP-conserving limit. 
One is that the additional mixing between the CP-odd and the CP-even components affects 
the significance of each process. The other is that Mf 23 are no longer independent of one 
another but constrained by (TT^]) . In combination with the latter relation, the mass order- 
ing ([H]) excludes the upper-right and the lower-left rectangular regions from the parameter 
domain. Another theoretical constraint that rules out the upper-left and the lower-right 
triangular regions, is perturbativity (fl~8]) . We have also checked numerically that imposing 
perturbativity guarantees the tree-level unitarity in Higgs-Higgs scattering in the parameter 
space under consideration. 

In addition to the above theoretical requirements, we impose the experimental constraints 
in ( )22|) in order to respect the current Higgs production limits from the LHC This rules 
out the dark-shadowed strip on the plane. The bounds on H 2 and H 3 apply as long as their 
masses are in the range between 160 and 600 GeV. All these conditions being fulfilled, the 
question to ask is how much portion of the remaining unshaded region the LHC can cover. 

As in the preceding plot, the W-H ± -Hi coupling largely determines the shape of the 
H ± Hi contours. As long as one moves along the horizontal axis, Hi is purely CP-even, 
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and therefore its interactions not involving other neutral Higgses reduce to those of h in the 
CP-conserving Higgs sectorJl For instance, the W-H ± -H\ coupling vanishes at the point 
(sin «i, sin a 2 ) = (sin/3,0) = (0.98,0), as is the case in figure l4al This explains why the 
significance decreases as one approaches this point. 

It should be meaningful to interpret the interesting behaviour of S around the point 
a i — a 2 — 0, although it belongs to the H 2 ^ exclusion strip. For instance, there might 
be a way to circumvent the constraint (|22l) by taking a different set of parameters. In the 
neighbourhood of this point, one finds a sudden drop of significance even though the W- 
H^-Hi coupling stays nearly maximal. This is because the Hi-b-b coupling is suppressed 
there as can be seen in table HI At the origin, H\ becomes fermiophobic, and its dominant 
decay modes are WW* and ZZ*. The fermiophobia of H\ generically gives rise to an 
enhanced B[H\ — > 77) ~ 1% without the charged Higgs loop contribution. Including the 
H ± -\oop, the diphoton branching fraction can range from zero up to (9(10%). [We give a 
related discussion in section |Vj] Therefore, the possibility of observing the H ± Hi state in 
the central low-significance hole via an alternative decay channel of Hi such as the 77 mode, 
is highly dependent on the other model parameters. 

As one departs away from the horizontal axis, sin«2 = 0, Hi acquires more of the CP- 
odd component. This reinforces the production rate of H ± H\ and thereby the significance, 
since the W-H ± -A coupling is not suppressed for any a± as can be seen in figure HH 

Overall, figures demonstrate that there is an ample portion of the parameter space in 
which the LHC has a good chance to discover the H ± Hi production. Obviously, the chance 
is substantially better in the "full luminosity" case, shown in figure l5bl With the exception 
of the narrow corners around both ends of the horizontal axis and the hole encircling the 
fermiophobic point, the discovery of H ± H\ through the rvbb final state should be possible 
at any point in the unshaded physically sensible domain defined by fl9]), (jTBl . and (122|) . 

In the same figures is presented the .Prpstrahlung process, which could also be detected 
depending on the mixing angles and the LHC operation condition. One can understand the 
shape of the W Hi contours by looking at the Hi-W + -W~ coupling in table [H In figure [Hal 
for "early searches", one finds no common portion of the parameter space with X > 5 for 
both processes. This setup predicts that it is not possible to discover both channels during 
the first data taking period in the Ivbb final states. Should this occur, this scenario will 
have to be modified to survive. In figure l5b~l one finds a much broader region in which 
ifi-strahlung could be discovered. Part of this region is shadowed in light grey where the 
"full LHC" could observe both this process and H ± H\ production. 



V. IMPACT OF iJ±-LOOP ON Hi 77 

This section is devoted to phenomenology of the lightest neutral Higgs. In particular, 
we show that the charged Higgs loop can exert a marked influence on the H\ — > 77 decay 
mode, which is of the utmost importance in the neutral Higgs search at the LHC. We shall 



2 Nevertheless, the interactions of H2 and H3 do not reduce to those of H and A due to non-vanishing 
sin «3. 
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(a) Pi (b) P 2 

FIG. 6. Cross-sections of the lightest neutral Higgs production mechanisms plotted against tan/3 
at y/s = 7 TeV for the benchmark points Pi and P%. The parameters are set as in (|20|) unless 
stated otherwise. The curves are independent of /i. 



consider the following standard Higgs production mechanisms: 

gg — > Hi (gluon-gluon fusion), 

qq — > jj + Hi (vector boson fusion), 

qq — >■ W/Z + Hi (Higgs-strahlung from W/Z), 

gg,qq — > ti+ Hi (associated production with a top-quark pair). 

In figures [61 we plot the cross-sections of these channels at y/s = 7 TeV against tan/3 for 
each of the two benchmark points. We have explicitly verified that all these production 
mechanisms are mostly unaffected in the range of Mi between 118 and 130 GeV and in the 
neighbourhood of the chosen M H ±. All the production rates are lower than those in the 
SM. This is because of the suppression of the couplings in table [B Combining production 
and decay, we obtain an intriguing result: Hi may have escaped from the past and present 
searches due to the suppression of both its production and the branching ratios to massive 
vector boson pairs. As this scenario generically predicts lower rates into fermion pairs, 



it could be disfavoured by the recent hint in the bb mode from the Tevatron [5JJ if it 
persists. Obviously, one can ease the reduction of Hi production by taking the Higgs mixing 
pattern to the SM-like limit, (sin a x , sin a 2 ) = (1,0). In this case, the charged-neutral Higgs 
associated production will be suppressed as shown in figures El Recall that we are assuming 
a hierarchical neutral Higgs mass spectrum and that all the experimental evidences from the 
LHC and the Tevatron are caused by the lightest of the three neutral Higgs particles, which 
is also supposed to be produced in association with H . Relaxing one or more of these 
hypotheses may lead to a different consequence, although we do not explore this possibility 
in this paper. 
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On the other hand, branching fractions of the two loop-induced decay modes, Hi — y 77 
and H\ — > gg, have possibilities of being enhanced as well as diminished, depending on the 
input parameters. If they are increased, the suppressed Higgs production can be alleviated. 

The decay into the digluon final state arises from the fermion loop, dominated by the 
top-quark contribution. Although the same type of loop is present also in the SM, its value 
can be different due to the variation of the Hi-t-t coupling which is determined by the 
Higgs mixing angles as in table HI In particular, an interesting role is played by the CP- 
odd component A in Hi that could enhance the digluon branching fraction. Its effects are 
twofold: (a) A does not couple to a vector boson pair unlike the CP-even Higgs h; (b) the 
digluon to difermion partial width ratio of A is higher than the corresponding quantity of 
h, i.e. 

^{A -)■ gg) = 2 3 T(h -» gg) 
T(A^ff) • T(h^ffY 

when Ma = Mh = 125 GeV. This effect stems from the different Lorentz structures of the 
scalar and the pseudo-scalar couplings with fermions. The potential enhancement of the gg 
decay mode might be of interest at a future lepton collider experiment. 

The decay amplitude of Hi — > 77 has two more components, the W- and the P^-loops, 
in addition to the fermion loop. In the SM, the dominant contribution comes from the 
W-loop. Like the fermion loop, the VF-loop can vary as a function of P^. 

One might attempt to increase the diphoton branching fraction by invoking the same 
mechanism as with the digluon mode. As Hi becomes pseudo-scalar-like, the fermion loop 
alone does receive the relative enhancement that is given by fT2"6]) with gg replaced by 77 and 
each numerator including only the fermion loop. However, this kills the jy-loop at the same 
time, which is the dominant contribution in the SM. Therefore, a high CP-odd fraction in 
Hi is not really helpful in increasing B(Hi — > 77). 

Nevertheless, it is known that the charged Higgs loop can greatly affect the 77 mode when 



M H ± is low |52H54| . In particular, this mechanism can lead to a high enough diphoton rate 
to be consistent with the excess surrounding the Higgs mass of 125 GeV recently observed by 
ATLAS [23[ and CMS j24[. The Pi-P + -P~ coupling depends on more parameters than the 
other loops, in particular on /1, with which one can play in order to engineer the amplitude. 
The maximal permitted size of the P^-contribution is mostly controlled by perturbativity 
and unitarity. In passing, we remark that change of /1 affects also r(P2,3 — > H + H~) and 
therefore the region excluded by (|22|) . as we will see in the following plots. 

To further study this potential variation, we plot the branching fraction of Hi — > 77 in 
figures [7a] and [TcJ Indeed, it can deviate substantially from its SM value depending on the 
mixing angles. One finds a valley between the two curves labelled "a" on each of the two 
plots, along which the branching ratio is highly suppressed. Across this valley, the Hi — > 77 
amplitude flips its sign. As the plots show, the contours are deformed if fi is changed, 
since they depend on the P ± -loop. One can also notice that part of the "d" contour of a 
high branching fraction lies on the verge of the perturbativity border, demonstrating that a 
sizeable P ± -loop contribution involves large quartic Higgs couplings. 

In conjunction with gluon-gluon fusion, this brings us to the cross-section of gg — > Hi — > 
77 in figures [7b] and [7dJ With our crude leading-order approximation, it is estimated to 
be about 20 fb in the SM. In the entire region of panel l7b| the cross-section is smaller 
than its SM value due to the suppressed Hi production. Nevertheless, one still has the 
option to change \i that affects the charged Higgs loop. For instance, the cross-section at 
(sin ai, sin 02) = (0.5,0) is approximately 6 fb for ji = 100 GeV. At this point, the W- 
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(a) Branching ratio, /i = 100 GeV 



(b) Cross-section, fi = 100 GeV 
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(c) Branching ratio, fi = 180 GeV 



(d) Cross-section, fi = 180 GeV 



FIG. 7. Contours of the branching ratios and the cross-sections involving the diphoton final state 
on the (sincvi, sin 02) plane at ^/s = 7 TeV. The two benchmark points Pi and P2 are marked by 
dots. The parameters are set as in (|20|) unless stated otherwise. The medium-shadowed regions 
(red online) are excluded by theoretical constraints. The dark-shadowed region (blue online) is 
ruled out by current limits on H2 3. 



17 



and the if -loops interfere destructively. Taking /i = 180 GeV instead, one can make the 
interference constructive by flipping the sign of the if ^contribution. This pushes the cross- 
section up above the SM prediction, as displayed in figure [7dl In this plot, the rate of 
gg — > Hi 77 can be comparable to its SM value even when the mixing angles are fairly 
removed from the SM-like point, (sina^, sina 2 ) = (1,0). In combination with figures 
this exhibits the exciting possibility that the LHC can discover a light charged Higgs in the 
foreseeable future while the recent evidences for Higgs decays to 77 are partly due to the 
charged Higgs loop effect. 

As one can expect, the charged Higgs loop effect weakens for high M H ±. For instance, 
B[H\ — > 77) for M H ± ~ 135 GeV is roughly half its value for M H ± = 86 GeV, when one 
keeps (sin «i, sin a 2 ) = (0.5,0) and /i = 180 GeV. 



VI. CONCLUSIONS 

We demonstrated the superb potential of the LHC to discover a light charged Higgs 
close to the jy-boson mass in the CP-violating 2HDMI. For this purpose, we exploited 
the property of the pp — > W — > H\H ± process that its cross-section is stable against 
variation of tan/3, which can be anywhere above the lower bound placed by the most stringent 
FCNC constraints. We set the Higgs masses to their favoured values, Mh± ~ 86 GeV and 
Mi ~ 125 GeV, and found that the LHC can spot an excess in the rvbb events, except in 
small corners and holes of the parameter space. We also examined the iix-strahlung channel 
from W, which offers access to the complementary parameter volume to that of HiH ± . 

With regard to CP- violation, we observed interesting consequences of the scalar-pseudo- 
scalar mixing. In particular, B(H\ — > gg) can be considerably affected by the composition of 
Hi. This result might be relevant to Higgs scrutiny at prospective linear collider experiments. 

Motivated by the recent evidences for Higgs production at the LHC, we studied cr(gg — > 
Hi — > 77). Due to the Higgs mixing matrix elements appearing at interaction vertices, 
the production of Hi is generically suppressed in comparison to the SM. Nonetheless, the 
if ± -loop amplitude can make a constructive interference with the standard H^-loop to raise 
BiHi — > 77). Depending on the parameter choice, this effect can be significant enough to 
increase the diphoton production rate from gluon-gluon fusion up to or even above the SM 
prediction. Alternatively, one can opt to enhance the Hi production itself by letting the 
Higgs mixing pattern approach the SM-like limit. In this case, the other processes with Hi 
going into the fermion and the massive vector boson pairs are also restored, the price to pay 
being suppression of the H X H ± rate. 

Related to this last point, it may be beneficial to study further possibilities of the light 
charged Higgs scenario in which one makes assumptions different from those underlying the 
present article. 
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Appendix A: Gauge-fixing Lagrangian of the 2HDM 



Assuming the 't Hooft-Feynman gauge, we choose a basis for the <J>j by expanding the 
Higgs doublets as 



-i( C(S w + - s p H + ) \ 

*1 I VT_+r)\+i(cpz—spr)3) I 
V2 J 



-i(spw + + cpH + ) 

V2 



(Al) 



where w (z) is the Goldstone boson of W ± (Z). We remark that the 773 pseudo-scalar field 
is orthogonal to the neutral Goldstone boson z. 

As for the Goldstone boson mass spectrum, it is possible to find a convenient way to write 
the mass matrix. Having each $12 the same group representation of the SM Higgs, following 
the notation of 55] , in the gauge-Goldstonqfl basis we find the following representation of 
the co-variant derivative: 

(9 \ 
g 
g 



V - Vl ' 2 



(A2) 



V -g x ) 



While the vector boson (and ghost) mass matrix is my = T>\(T>\ t 
Goldstones mass matrix is: 

m^ = (I>i) T D 1 + (I> 2 ) T P 2 , 

therefore we get 




V 2 (V 2 ) T , the 
(A3) 



(A4) 



The mass matrix in Equation flA4[) shows that the Goldstones have a mass that is equiv- 
alent to the SM-ones, as expected since the gauge sector has not been extended. 

As we have already intimated, the ghost mass matrix and interactions are defined by 
means of the same matrix T> via 



rn 



ghost 



V 1 (V 1 ) T + V 2 (V 2 ) 



(A5) 



Notice that the m 2 ghost and the m 2 , of equation (IA4I) have different numbers of zero- 
eigenvalues, but their non-zero eigenvalues are in a one-to-one correspondence; furthermore, 
the eigenvalues of the gauge- fixing mass matrix are the same of the gauge boson mass matrix, 
as expected. 

Then, the ghost Lagrangian is defined as 



C 



ghost 



C" 



V a 2 ■ (V\ + S b 2 ) 



T 



(A6) 



3 The 4x3 matrix follows from the four gauge bosons Z, and the three Goldstone bosons 4> l \i = i^. 
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where the matrices Si >2 represent the link between the fluctuations (Goldstones) of the 
Higgses around their VEVs and the gauge bosons; a convenient^] way to write this matrices 



is 



where 



gr)i cpgz -c p gw 2 -cpgiw 2 

(SiY = 7; ( -cpgz gm c p9 w i cpgiwi | > ( A7 ) 

cpgw 2 -cpgw x gr} x -g^ 

gr] 2 spgz -spgw 2 -spg x w 2 
{S 2 ) T = ^ [ -spgz gr] 2 s f} gw 1 spg^ | , (A8) 
spgw 2 -spgwY gr] 2 -gity 



w + + w ... 
wi = 7=—, (A9) 



V2 

7! 



w 2 = i -= — . (A10) 



Finally, the ghost fields ( ( c) read as 



c= (wi w 9 2 B g ) . ( A11 ) 

where 

W f = -JlJ_JL, (A12) 
— w~ 

< = ( A] - 3 ) 

wf = cos^iy% + sm.O w A g , (A14) 

fi 9 = -sin^vy-Zg + cosQ w A g . (A15) 



Appendix B: Perturbativity at high tan/3 

In order to properly realise the high tan /3 scenario, we must carefully consider the per- 
turbativity of the couplings. Firstly, from Equation (1T2"|) we see that the tan/3 — > 00 ap- 
proximation leads to the following definition of the heaviest neutral scalar mass: 

From Equations ffTTj) we can conclude that all the couplings are finite and perturbatively 
small except for Ai and A3JI In fact, for tan /3 — > 00 or equivalently cos /3 — > and sin /3 — > 1, 



4 In fact, the explicit calculation of each single matrix gives rise to connection between ghosts and both the 
charged Higgses and the 773 combination of scalar fields. Consistently, the sum of the two contributions 



disarms these connections: ghost fields do not link to charged Higgs. 
5 While also ImAs is apparently divergent, by plugging the definition of Af| (tan j3 — > 00) in (|17f|) it is easy 
to see that the result is always finite. If v\ exactly vanishes, one can rotate away the phase of A5 to render 



CP-conservation manifest 56 
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FIG. 8. Differential cross-sections of the process pp — > W — >■ H H\ — > r^ubb with respect to 
of the 6-quark and the r-lepton at the indicated centre-of-mass energies, for the benchmark point 
P\. Each curve is the sum of da(pp — > r + vbb)/dpT and that for the conjugate final state. The 
shadowed range (red online) shows the selection cuts set for the signal-to-background analysis. 



the two potentially diverging couplings become: 

1 



Ai 



A, 



v 2 cos 2 j3 
Ml - M 2 2 R 22 R 12 
v 2 cos /3 -R33 



-K32-K33 



(B2) 
(B3) 



Hence, no matter how large tan/3 is, we can always keep Ai and A3 perturbatively small by 
choosing 



-R32-R33 



Mi ~ M 2 , 



(B4) 
(B5) 



respectively. This shows that there is no strict upper limit on tan/3 from perturbativity. 
However, we do not make use of this mechanism for drawing our conclusions in the main 
part of the present paper. 



Appendix C: pt distributions of b and r 



In figure El we show the pr spectra for the 6-quark and the r-lepton, extracted from the 
process pp — > W — > H Hi —> r^isbb at the LHC. They were evaluated at the benchmark 
point Pi. In the px,b distributions, one finds that the selection cuts in fl25l) act on regions 
with suppressed differential cross-sections. This means that we retain the major portion 
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of events at our disposal even after applying the cuts. Even in the pr, T distributions, the 
excluded area below each curve is relatively small compared to the total cross-section. 
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